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The inverse mass expansion of the one-loop effective action of a scalar particle 
with mass m in the background of a (non Abelian) gauge field A^ and a matrix 
valued scalar potential V is usually obtained from the one-loop determinant of 
the corresponding fluctuation operator M, which (in Euclidean space) reads 

M = -D"^ + m^ + V{x) with D^ = d^, + igA^{x) . (1) 

by using the Schwinger proper time representation, 

TeffiA, V] = -ln(detM) = -Tr(lnM) = / — Tr e"™ , (2) 

JO J 

and expanding in powers of the proper time parameter T. For the gauge case 

the coefficients of this expansion were calculated up to 0(T^) by field theory 

methods in [1]. The operator trace is nothing but the diagonal element of 

the heat kernel for the operator M, integrated over space-time. Using standard 

(recursive and non-recursive) heat kernel methods [2] the coefficients have been 

computed up to 0{T ) for the pure scalar case and 0{T^) for the gauge case 

only recently [3]. 

The worldline approach to the inverse mass expansion makes use of progress in 

calculations of one-loop amplitudes by string inspired techniques [4, 5, 6]. 

The operator trace is written as a one-dimensional path integral over the space 

of closed loops in space-time with fixed proper time circumference T. The 

one-loop effective action takes the form 



T[A,V]= ^e-™^trP/ Vxexp - dri^ + igx'^A^ + V) . (3) 

Jo 1 Jx(T)=x(0) L Jo 4 J 



Here V denotes path ordering and the operator trace has reduced to an ordinary 
matrix trace. After expansion of the interaction exponential the path integral 
can be evaluated by Wick contractions. The appropriate Green function is the 
Green function of the Laplacian on the circle with periodic boundary condi- 
tions [5]. Since the defining equation of this Green function has no solution, 
one is forced to introduce a 'background charge' p on the worldline. A very 
convenient choice is a uniformly distributed background charge p = 1/T (which 
distinguishes the present formalism from former path integral approaches in 
this context [7]), yielding the translationally invariant Green function 

GB{Tl,T2) = \n-T2\-^^^^P^. (4) 

Due to the existence of the zero mode this Green function cannot be applied to 
the path integral as it stands. One has to introduce a loop center of mass xq 
and a relative coordinate y, 

x'^(r) = x(; + y^ir) with / dr y^ir) = , (5) 

Jo 

and separate the integration over the center of mass from the path integral. 



jVx= I d'^xo I Vy 



(6) 



The resulting path integral is over the space of all loops with a common center 
of mass xq. Using y^ = x^ this yields 



T[A,V] 



dT 



-ra^T. 



T y2 



tr V /d% /2?yexp - / dT{^+igr A^+V) . (7) 



If one chooses Fock-Schwinger gauge for the background gauge field with refer- 
ence to the center of mass xq [6] , it can be written as 



(8) 



A^{xQ + y) = yP dr] T] Fp^ixo + r]y) 
Jo 

and the background fields F and V can be expanded covariantly: 



Fp^{xo + vy) = e'>''''F,^{xo), V{xo + y) = e^^Vixo) . (9) 

Using these formulas the expansion of the interaction exponentials in Eq. (3) 
in the proper time parameter T yields the following manifestly covariant form 
of the action: 



nF,V] = r^e-'^trfd'^xof:^~^ 







T 



n=0 



n 



Vy exp 
1 



\X 



T1=T rTn-1 

dr2 . . . / dTn 
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[xo, 



.(10) 



The factor - arises upon eliminating of the first r-integration, which is possible 
due to the freedom of choosing the origin somewhere on the loop. 
For the calculation of the heat kernel coefficients to a given order in T one has 
to perform the following steps: 

• Wick contractions 

The sum in Eq. (10) is truncated and the exponentials are expanded. Then 
all possible Wick contractions have to be evaluated using the contraction 
rules 



(2/^(ri)y'^(r2)) 
{y''iri)y''ir2)) 
{rin)rir2)) 



-9^''GBin,T2) = -g 



fiu 



Tl -T2 



in - T2) 



-9^''Gb{ti,T2) = -gf"" sign(Ti - T2 



T 

2(ti - T2) 



+9^''GB{n,T2) = +g'"'[25iTi -T2)-- 
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where the dot denotes differentiation with respect to the first variable. 
From these rules one may alternatively derive contraction rules for expo- 
nentials, which can be used equivalently. 



• Integrations 

The polynomial r- and r/-integrations are performed. After the expansion 
in the first step the r/-integrations are trivial. The integrands of the r- 
integration involve the worldline Green function Gb and its derivatives. 
They are polynomial in the variables Tj. 

• Basis reduction 

The number of terms will be drastically reduced by an appropriate basis 
reduction, which is essential for the practical use of the coefficients in 
numerical calculations. In the pure scalar case there is a unique basis 
with the property that no box operators □ = 5^ occur. With the choice 
of the Green function Eq. (4) no self contractions exist, i.e. box operators 
do never occur and no partial integrations have to be performed to reduce 
the result. The only operation left is the identification of cyclic equivalent 
terms, which is trivial. The gauge case is much more complicated and 
requires a reduction algorithm, which is described below. 

Finally one obtains the result in the form 

nF,V] = r^[4nT]-'^/'e--'^f:tlLfdxotrOn . (12) 

Jo T ^1 "! -^ 

The factor [47rT]^'^'^ arises from the normalization of the free path integral. 
For the scalar case the calculation of the coefficients was done up to On [8, 9]. 
In the following we describe the basis reduction algorithm for the gauge case as 
proposed by Miiller [10]. It turns out that there is a minimal basis of invariants 
without box operators and partial integrations are still not necessary to reduce 
our results into this basis. However, besides the cyclic permutations one has 
to use the antisymmetry of the field strength tensor, the Bianchi identities and 
the exchange of covariant derivatives, 

D^D^X = D^D^X + ig [F^^,X\ , etc. (13) 

for a further reduction. 

This is done by the following steps (for a more detailed discussion of the algo- 
rithm and a proof of minimality see [10]): 

• Elimination of 'middle' derivatives 

The contractions of derivatives of a field strength tensor belong to different 
classes with respect to the contractions of the field strength tensor itself. 
This can be seen very easily in a diagrammatic picture, where the (cyclic) 
function tr is represented by a circle, the fields by points on the circle and 
the contractions between the fields by lines crossing the circle. In this 
picture the indices of a field strength tensor divide the circle into three 
sectors (left, middle, right). Consequently the derivatives fall into different 
classes, according to the sector, they are contracted with. The Bianchi 
identity involves all classes of derivatives and can be used to eliminate 
one of them. It is useful to take the symmetric choice and eliminate the 



middle derivatives. This has to be done with all derivatives of a field 
strength tensor, which possibly requires exchange of derivatives. 

• Reduction of multiple contractions between factors 

Terms with multiple contractions are reduced by the following rules 



tr^. . . r ^1/ . . . L/^t'iiK. 



tr(. . . D^F^,^ . . . D^F^ 



^x 



tT{...F^,...D^D,X...) 



^tv{...F^,...D^F,^...), (14) 

-tT{...D^F^^...D^F^^...) 

+ tTi...D^F,^...D^F,x...), (15) 

ii<7tr(...F^....[F^.,X]), (16) 



which involve the Bianchi identity, the antisymmetry of the field strength 
tensor and the exchange of derivatives. The aim of these rules is to pro- 
duce contractions of the field strength tensors among themselves whenever 
possible. 

• Final arrangement of indices and factors 

The indices still have to be ordered using the antisymmetry of the field 
strength tensor and the exchange of derivatives using e.g. a minimal con- 
traction distance rule. After this the cyclic equivalent invariants are re- 
duced as in the pure scalar case. 

The reduction algorithm produces invariants with smaller number of deriva- 
tives and higher number of field strength tensors by exchange of derivatives. 
Therefore one has to start the reduction with the terms containing a maximum 
number of derivatives. 

For a more specialized situation one can use existing additional symmetries, like 
the mirror symmetry [10], for a further reduction. For the general case Table 
1 gives an overview of the number of invariants in the minimal basis to given 
order and number v of occurring scalar background fields. 
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Table 1: Number of basis invariants in different orders of the expansion. 

The computerization of the calculation splits into two parts. The first consists 
of the expansion of the interaction exponentials and the Wick contractions. As 



in the scalar case this is conveniently done with FORM [11] up to 0{T^). The 
second part of the calculation consists of performing the r-integrals and the 
application of the basis algorithm to the raw coefficients. The coding of the 
basis algorithm requires a rule based system for symbolic manipulation with 
fast and flexible pattern matching. For this purpose we chose the new language 
'M' developed by P. Overmann [12]. The reduction has been completed for the 
coefficients up to O(T^), the reduction of Oq is currently done. Up to 0{T^) 
the results have been checked to be equivalent with the result obtained from a 
modified non-recursive heat kernel method [13]. Additionally, a check with the 
results of [1] has been done up to 0(T^). 

The expression for O5 is too large to be presented here and can be found in [9]. 
The results to 0{T'^) read (absorbing the coupling constant g into the fields, 
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In conclusion, we have obtained the inverse mass expansion of the one-loop 
effective action for a massive scalar particle in the background of a non Abelian 
gauge field and a matrix valued scalar potential up to 0{T^), reduced to a 
minimal set of unique basis invariants. Within the framework of the worldline 
path integral method it is currently feasible to compute the coefficient Og . Since 
the basis reduction of this coefficient has not been completed yet, the results 
will be presented in a forthcoming publication [13]. 
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